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ABSTRACT 


An  experiment  was  conducted  which  measured  the  Radio  Fre¬ 
quency  (RF)  radiation  from  the  PHERMEX  accelerator,  capable  of 
30  MeV  and  600  A.  This  was  accomplished  by  placing  TEM  horn 
antennae  at  varying  angles  from  the  path  of  the  electron  beam.  The 
signals  received  by  the  antennae  were  then  recorded  by  using  a  Digi¬ 
tizing  Camera  System  (DCS).  Measurements  were  taken  of  the  radia¬ 
tion  from  propagating  and  non-propagating  beams,  beams  with  energy 
above  and  below  Cherenkov  threshold,  and  beams  with  varied  cur¬ 
rents.  The  captured  RF  signals  and  their  corresponding  frequency 
spectra  were  then  analyzed.  This  analysis  showed  that  the  radio  fre¬ 
quency  radiation  from  the  beams  below  the  Cherenkov  threshold 
contained  primarily  transition  radiation;  when  above,  diffracted 
Cherenkov  radiation  was  observed.  Non-propagating  beams  produced 
larger-angle  radiation  and  had  less  definition  in  their  spectrum.  All 
electric  fields  measured  were  proportional  to  the  beam  current. 
Lastly,  the  electron  beam  pulse  width  and  separation  were  determined 
by  both  the  received  signals  and  their  spectrum. 
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I.  INTRODUCTION 

A.  BACKGROUND 

1.  Historical  Look  at  Cherenkov  Radiation 

The  following  history  of  Cherenkov  radiation  discovery  and 
development  is  based  on  Jelley's  introduction  [Ref.  1]  and  follows  the 
summary  by  O’Grady  [Ref.  2]. 

As  early  as  1910,  the  effects  of  Cherenkov  radiation  were 
observed  by  Mme.  Curie.  Her  bottles  of  concentrated  radium  solution 
emitted  a  pale  blue  light.  This  phenomenon  was  subsequently  partially 
investigated  by  Mallet  during  the  years  between  1926  and  1929.  How¬ 
ever.  it  was  not  until  1938  that  Cherenkov  commenced  an  exhaustive 
series  of  experiments  to  fully  understand  this  phenomenon.  During  this 
time,  a  theory  was  proposed  by  Frank  and  Tamm  which  had  excellent 
agreement  with  Cherenkov’s  experimental  results.  Since  their  work,  there 
have  been  numerous  theories  and  experiments  concerning  Cherenkov 
radiation  detection  and  Cherenkov  radiation  applications. 

By  Frank  and  Tamm’s  definition  [Ref.  3],  Cherenkov  radiation 
takes  place  when  a  charged  particle  moves  at  a  constant  velocity  which  is 
larger  than  the  phase  velocity  of  light  in  a  transparent  medium.  Electrons 
can  be  accelerated  to  speeds  greater  than  the  velocity  of  light  in  air  and 
thus  produce  Cherenkov  radiation.  Frank  and  Tamm’s  theory,  however, 
contains  a  number  of  simplifying  assumptions.  Two  of  the  more  impor¬ 
tant  assumptions  are  that  the  medium  is  an  unbounded  continuum  and 


1 


the  track  length  is  infinite  [Ref.  1].  Applying  these  assumptions  to  most 
accelerator  applications  is  unrealistic.  When  an  electron  bunch  leaves 
the  vacuum  of  an  accelerator,  another  medium  (such  as  air)  is  encoun¬ 
tered.  The  electron  beam  steadily  loses  energy  due  to  ionization  of  air 
molecules.  Therefore,  the  beam  has  a  finite  length  and  propagates  in  a 
bounded  medium. 

When  the  two  previously  made  assumptions  are  violated,  two 
similar  types  of  radiation  occur— diffracted  Cherenkov  radiation  and 
transition  radiation  [Ref.  4].  Diffracted  Cherenkov  radiation  is  caused  by 
the  finite  beam  interaction  length  and  transition  radiation  is  caused  by  a 
boundary  where  the  index  of  refraction  of  the  medium  changes.  Another 
type  of  radiation  closely  related  to  transition  radiation  is  diffracted 
transition  radiation  [Ref.  51.  In  this  case,  radiation  is  produced  when  a 
charged  particle  passes  through  an  aperture  or  travels  near  any  interface 
between  differing  media.  The  circular  beam  port  of  an  accelerator  is  an 
obvious  source  of  diffracted  transition  radiation.  However,  since  the 
circular  beam  port  dimensions  are  constant,  the  resulting  radiation  will 
merely  be  added  to  the  transition  radiation  and  considered  the  electro¬ 
magnetic  pulse  (EMP). 

Since  the  major  effects  of  this  experiment  concern  both  EMP 
and  diffracted  Cherenkov,  a  more-detailed  summary  for  each  is  contained 
in  the  following  two  sections.  The  last  section  in  this  chapter  addresses 
the  similarities  and  differences  of  both. 
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2.  Cherenkov  and  Diffracted  Cherenkov  Radiation  Theory 

As  previously  stated,  Cherenkov  radiation  is  produced  whenever 
a  charged  particle  travels  through  a  dielectric  medium  faster  than  the 
phase  velocity  of  light  in  the  same  medium.  The  underlying  reason  radia¬ 
tion  forms  is  that  the  charged  particle  temporarily  excites  the  medium. 
The  ionized  molecules  briefly  radiate  as  a  dipole.  Of  particular  interest  is 
the  Cherenkov  angle  0C.  When  the  interaction  length  is  infinite,  the  radia¬ 
tion  propagates  only  at  this  angle.  The  Cherenkov  angle  can  be  calcu¬ 
lated  using 

cos  0C  (1) 

where  c  is  the  speed  of  light  in  a  vacuum,  n  is  the  refractive  index  of  the 
medium,  and  v  is  the  speed  of  the  particle.  A  more-detailed  description  of 
Cherenkov  radiation  and  a  thorough  derivative  of  the  fundamental  equa¬ 
tion  for  output  of  radiation  is  given  by  Jelley  [Ref.  1J.  However,  since  0C 
occurs  at  approximately  one  degree,  the  diffracted  Cherenkov  radiation 
will  be  more  applicable  to  this  experiment. 

Measurements  by  Neighbours,  Buskirk,  and  Saglam  [Ref.  6]  in  1983 
showed  that  the  Cherenkov  cone  angle  is  broader  as  the  interaction 
length  is  decreased  from  infinity.  This,  of  course,  is  diffracted  Cherenkov 
radiation.  Neighbours,  Buskirk,  and  Saglam  [Ref.  7]  derive  an  equation 
for  the  total  coherent  power  per  unit  solid  angle,  radiated  at  the  fre¬ 
quency  v  by  a  periodic  charged  particle  beam  traveling  a  finite  distance  L 
at  constant  velocity.  This  equation  is 
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w(v,k)  =  v*QR* 


(2) 


where  v0  is  the  fundamental  frequency  of  the  beam  generator  and  where 
the  constant  Q  is 


where  \x  is  the  permeability  of  the  medium  and  q  is  the  change  in  the 
electron  bunch.  The  radiation  function  R  is  given  by: 

R  =  witn  sinQ  I(u)F(k)  (3) 

where: 

6  -  angle  between  emitted  radiation  and  beam  path 

T)  =  LA 

L  =  beam  Interaction  length 

X  =  radiation  wavelength 

F(k)  =  dimensionless  form  factor 

I(u)  =  sin  u/u  =  diffraction  function 
where  u  is: 


u  =  Tcqfcos  0e  -  cos  0]  (4) 

A  graph  of  equation  3  with  the  appropriate  constants,  for  this  experiment 
yields  Figure  1.  This  envelope  function  determines  the  relative  output 
power  for  a  given  angle  from  the  beam  path. 
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Figure  1.  Envelope  Function  for  Beam  Parameters 

Two  additional  comments  concerning  diffracted  Cherenkov  radi¬ 
ation  are:  (1)  “the  power  emitted  for  coherent  radiation  is  proportional  to 
the  source  of  the  beam  current”  [Ref.  81.  (2)  Neighbours,  Buskirk,  and 
Maruyarna  [Ref.  9]  also  predict  that  “the  oscillations  of  the  frequency 
distribution  of  the  radiation  are  of  increasing  "frequency”  as  the  angle  of 
propagation  moves  toward  the  back  direction.”  To  better  understand  this 
statement,  one  should  view  Figure  2  (obtained  from  Reference  9).  From 
Figure  2,  the  predicted  pattern  is  clear. 
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Figure  2.  Frequency  Shifting  with  Angle 

3.  Transition  Radiation  Theory 

While  Cherenkov  radiation  occurs  when  the  velocity  of  the  par¬ 
ticle  is  greater  than  the  speed  of  light  in  that  medium,  transition  radia¬ 
tion  can  occur  whether  the  beam  energy  is  above  or  below  the  Cherenkov 
threshold  energy.  When  the  electron  passes  through  a  boundary  between 
two  different  media  transition  radiation  is  produced.  Frank  and  Ginsburg 
[Ref.  10]  explain  the  nature  of  this  phenomenon  by  considering  the 
transition  of  an  electron  from  a  vacuum  into  a  perfect  conductor.  While 
the  electron  moves  in  the  vacuum,  there  is  no  radiation.  However,  the 
field  in  the  vacuum  is  equal  to  the  field  of  the  electron  and  its  image 
moving  toward  it.  Once  the  electron  exits  the  vacuum,  both  the  electron 
and  its  image  cease  to  exist  from  the  vacuum  field  perspective.  Therefore, 
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radiation  must  be  given  off  in  much  the  same  way  as  it  is  emitted  when 
an  electron  suddenly  stops. 

Work  done  by  Busldrk  and  Neighbours  [Ref.  4]  develops  a  useful 
equation  by  which  the  electric  and  magnetic  field  strengths  can  be  evalu¬ 
ated.  The  maximum  magnetic  field  is  given  by 

B  =  - 1 - 

H  R  l-pcos0 

where  0  is  the  angle  between  the  beam  path  and  detection  point,  R  is  the 
distance  between  the  field  point  and  the  beam  port,  and  I0  is  the  current 
of  the  accelerator.  The  constant  (3  is  defined  to  be 


where  c  is  the  velocity  of  radiation  in  the  median,  cD  is  the  velocity  of 
light  in  free  space,  and  n  is  the  index  of  refraction  of  the  medium.  For 
this  experiment,  the  transition  will  be  from  vacuum  to  air  as  the  electron 
bunch  leaves  the  accelerator.  The  electric  field  can  therefore  be  written  as 
[Ref.  11] 


E  =  £S2i2Mi _ 1 _  (6) 

4tc  R  °1-Pcos0 

where  co0  is  the  permeability  of  free  space.  Figure  3  [Ref.  11]  shows  the 
relation  between  transition  energy  emitted  per  angle  as  a  function  of 
beam  energy. 
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Figure  3.  Transition  Energy  Emitted  per  Unit  Angle 

Diffraction  radiation  during  this  experiment  will  be  measured 
and  treated  as  additional  transition  radiation. 

4.  Comparison  of  Diffracted  Cherenkov  and  Transition 

Radiation 

The  similarities  between  Cherenkov  and  transition  radiation 
were  summarized  by  Buskirk,  Mack,  and  Neighbours  [Ref.  12].  First, 
both  radiation  fields  are  polarized,  with  the  electric  field  in  the  plane  of 
the  observer  and  the  beam  path,  and  the  magnetic  field  perpendicular  to 
that  plane.  Also,  the  electric  field  pulse  qualitatively  follows  the  current 
pulse.  Last,  the  field  from  either  diffracted  Cherenkov  or  transition  is 
proportional  to  the  beam  current.  Therefore,  the  radiated  power  is  pro¬ 
portional  to  the  current  squared. 
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There  are,  however,  some  important  differences  between 
Cherenkov  and  transition  radiation.  Transition  radiation  is  concerned 
with  transiting  boundaries.  Boundaries  associated  with  electron 
accelerators  are  the  beam  port  and  beam  stops.  Although  data  was  taken 
with  beam  stops  in  position,  they  will  be  analyzed  in  a  follow-on 
experiment.  The  next  difference  is  that  Cherenkov  radiation  occurs  in  air 
when  the  beam  energy  is  above  the  Cherenkov  threshold.  Also,  the 
Cherenkov  power  radiated  (without  diffraction)  is  proportional  to  the  path 
length.  Figure  4  [Ref.  4]  shows  two  regions:  A  and  B.  Cherenkov  radiation 
(without  diffraction)  is  concentrated  at  one  angle  (region  B).  With  diffrac¬ 
tion.  radiation  spreads  out  about  the  Cherenkov  angle  0C.  whereas  tran¬ 
sition  radiation  occurs  at  all  angles  (regions  A  and  B)  with  increased 
intensity  in  the  forward  direction. 

B.  PREVIOUS  EXPERIMENTS 

There  have  been  a  number  of  experiments  conducted  which  measure 
optical,  x-ray.  or  microwave  radiation  from  electron  beams.  However,  rel¬ 
atively  few  experiments  have  measured  radio  frequency  radiation  from 
propagating  electron  beams.  Buskirk  and  Neighbours  have  previously 
performed  a  number  of  such  experiments  [Ref.  13)  dealing  with  radio 
frequency  radiation.  The  present  experiment  differs  from  their  earlier 
efforts  in  the  following  ways:  (1)  the  number  of  angles  at  which  data  were 
sampled  was  Increased,  (2)  many  measurements  were  recorded  with 
varied  beam  parameters,  and  (3)  most  importantly,  a  digitizing  camera 
was  incorporated  into  the  collection  system.  This  provided  a  quicker  and 


9 


easier  method  for  data  collection  and  storage.  This  digitizing  camera  also 
has  another  important  advantage  over  the  previous  analog  system.  The 
resulting  data  files  can  be  analyzed  by  computer  programs  and  software 
without  painful  conversions. 


C.  PURPOSE 

The  general  overall  objective  is  similar  to  that  of  earlier  experiments: 
to  understand  the  radio  frequency  radiation  produced  by  an  intense  elec¬ 
tron  beam.  More  specifically,  the  purpose  of  this  experiment  can  be 
defined  in  the  following  statements: 

1.  Observe  differences  in  propagating  to  non-propagating  beams  by 
varying  the  beam  focus. 

2.  Attempt  to  isolate  diffracted  Cherenkov  radiation  effect  from  elec¬ 
tromagnetic  pulse  (EMP)  by  varying  the  beam  energy  above  and 
below  the  Cherenkov  threshold  (about  20  MeV). 

3.  Observe  possible  instabilities  in  beam  propagation  due  to  high  cur¬ 
rents.  Conversely,  investigate  the  effect  low  current  has  on  pinching 
the  beam. 

4.  Analyze  the  radiation  signal  and  frequency  spectrum  for  empirical 
observations. 
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n.  EXPERIMENTAL  APPROACH 


A.  EQUIPMENT  ARRANGEMENT 

The  specific  equipment  used  in  this  experiment  is  listed  in  Appendix 

A.  In  order  to  collect  data  simultaneously  from  two  antennae  and  also 
obtain  the  original  source  signal,  three  independent  collection  systems 
were  required.  The  equipment  arrangement  for  one  of  the  antennae  col¬ 
lection  systems  can  be  seen  in  Figure  5.  The  collection  system  for  the 
source  signal  is  essentially  the  same,  except  that  one  of  the  accelerator’s 
B  dot  loops  at  the  beam  port  replaces  the  antenna.  Actual  pictures  of  the 
collection  systems  (minus  the  antennae)  are  seen  in  Figure  6.  In  this 
setup,  the  radiation  from  the  electron  beam  is  received  by  the  TEM  horn 

*  antenna  and  transmitted  to  the  oscilloscope  by  means  of  the  doubJe- 
shielded  transmission  cable.  The  signal  on  the  oscilloscope  is  digitized  by 

* 

a  Digitizing  Camera  System  (DCS)  and  stored  on  a  computer.  The  B  dot 
loop  signal  was  recorded  in  a  similar  fashion  and  was  also  used  to  trigger 
the  other  two  oscilloscopes. 

B.  ACCELERATOR  AND  BEAM  CHARACTERISTICS 

The  accelerator  used  to  produce  the  electron  beam  was  PHERMEX, 
located  at  Los  Alamos  National  Laboratory  (LAND,  New  Mexico.  A 
schematic  of  PHERMEX  (Figure  7)  shows  the  three  cavities  which  are 
excited  by  means  of  radio  frequency  amplifiers  operating  at  50  MHz. 

* 

PHERMEX  can  produce  a  30  MeV,  900  A  beam  focused  to  a  6  mm 
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Figure  7.  Schematic  of  Phermex  Accelerator 


diameter.  The  energy,  current,  and  diameter  can  be  varied  as  required 
(see  Beam  Adjustment  section).  The  pulse  train  created  by  PHERMEX  is 
a  maximum  of  10  pulses,  two  to  three  nanoseconds  (ns)  long,  and  sepa¬ 
rated  by  20  ns.  Pul3e  trains  can  be  formed  up  to  one  per  minute  [Ref. 
14). 


C.  ANTENNA  CHARACTERISTICS 

The  transverse  electromagnetic  (TEM)  horn  antenna  was  developed 
at  the  National  Bureau  of  Standards[Ref.  15).  The  horn  is  configured  as  a 
two-conductor,  end-fire,  traveling-wave  structure.  Continuously  tapered 
resistive  loading  of  the  two  conductive  elements  provides  a  relatively 
short,  broadband,  nondispersive  antenna  which  has  high  directivity. 
Sections  of  polystyrene  foam  support  the  horn's  reslstivcly  tapered  con¬ 
ductors.  The  antenna  is  enclosed  in  a  bakelite  box  which  measures 
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about  30  cm  wide  by  15  cm  high  by  51  cm  long.  A  schematic  of  the 
antenna  is  shown  in  Figure  8.  The  antenna  can  measure  electric  fields 
between  1  and  1.600  from  1  MHz  to  15  GHz,  with  a  flatness  of  ±2  dB 
(Ref.  16]. 


Figure  8.  Schematic  of  TEM  Horn  Antenna 


D.  DESCRIPTION  OF  DIGITIZING  CAMERA  SYSTEM  (DCS) 

The  Tektronix  DCS  (Ref.  17]  is  a  charge-coupled  device  (CCD)  cam¬ 
era  which  connects  directly  to  the  face  of  the  oscilloscope.  490  vertical  by 
384  horizontal  picture  elements  are  used  to  form  an  image  of  the  trace  on 
the  oscilloscope.  Then  the  camera  transfers  the  analog  video  signal  to  the 
frame  store  board.  The  frame  store  board  is  a  General  Purpose  Informa¬ 
tion  Bus  (GPIB)  located  in  an  IBM-compatible  computer.  The  conversion 
from  analog  video  signal  to  digital  format  is  performed  by  the  frame  stor¬ 
age  board.  Here,  the  digitized  raw  image  which  is  held  in  the  computer’s 
temporary  memory  can  then  be  stored  on  diskettes  or  a  hard  disk.  The 
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DCS  software,  along  with  the  frame  store  board,  allows  for  format  con¬ 
version.  wave  form  manipulation,  and  file  transfer.  Williams  provides  a 
more-detailed  summary  of  the  DCS  system  and  software  [Ref.  181. 

E.  ANTENNA  PLACEMENT 

The  two  antennae  were  each  mounted  on  top  of  eight-foot  ladders.  At 
this  height,  the  antennae  were  approximately  two  and  a  half  feet  higher 
than  the  beam  port.  Figure  9  shows  one  of  the  actual  antenna  place¬ 
ments.  However,  in  this  picture,  both  antennae  are  mounted  on  the  same 
ladder.  The  positioning  of  each  antenna  was  relative  to  the  beam  path 
from  the  beam  port.  At  a  constant  radius  of  90  feet,  six  positions  (5-,  10-, 
15-,  30-,  45-,  and  60-degree  angles)  relative  to  the  beam  path  were  mea¬ 
sured  and  marked.  A  schematic  of  the  antenna  placements  is  seen  in 
Figure  10.  Permanent  structures  on  one  side  of  the  beam  path  prevented 
measurements  in  that  direction. 

F.  DATA  AC9UISITION 

1.  Beam  Adjustment 

Radiation  from  the  electron  beam  was  recorded  in  sets  of  three 
to  five  signals  per  antenna  angular  position.  This  procedure  was  carried 
out  to  form  specific  groups  comprised  of  four  variations  of  the  electron 
beam  (see  Table  1).  The  first  group  of  data  sets  was  taken  at  full  current 
(510  A)  and  full  energy  (29.5  MeV).  The  beam  was  focused  to  6  mm  for  all 
shots  in  this  group.  The  second  group  consisted  of  data  sets  received 
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Baam  Path 


Figure  10.  Antenna  Placement 

TABLE  1 

DATA  COLLECTION  SUMMARY 


Angles  (Degrees) 


Group 

Beam  Parameters 

Antenna 

5 

10 

15 

■71 

45 

T  •  ^  ' 

1 

Full  Current  & 
Full  Energy  & 
Tight 

A 

B1 

X 

X 

X 

X 

X 

X 

X 

X 

2 

Half  Current  & 
Full  Energy  & 
Tight 

A 

B1 

X 

X 

X 

X 

X 

X 

3 

Full  Current  & 
Full  Energy  & 
Diffuse 

A 

B2 

X 

X 

X 

H 

X 

X 

4 

Half  Current  & 
Half  Energy  & 
Tight 

A 

B2 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

with  the  beam  at  half  current  (300  A).  All  other  parameters  of  the  beam 
remained  fixed.  The  third  group  of  data  sets  was  altered  only  by  the 
diameter  of  the  beam.  The  accelerator's  focusing  magnet  was  adjusted  to 
deliver  a  diffused  or  under-focused  beam.  In  the  last  group,  a  half-energy 
(15  MeV),  half-current,  focused  beam  produced  the  signals. 

Throughout  all  of  the  shots  in  this  experiment,  the  pulse  train 
remained  at  either  six  or  seven  pulses  with  a  pulse  width  of  3.5  ns  and  a 
pulse  separation  of  20  ns  (see  Figure  11).  Therefore,  the  pulse  trains  var¬ 
ied  from  121  to  145  ns. 


SIGNAL  (Volts) 

Figure  11.  Typical  B  Dot  Signal 
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The  current  of  the  beam  can  be  determined  from  the  maximum 
peak  value.  The  beam  current  I  is  related  to  the  peak  voltage  V,  by 

I  *  CiV  (7) 

where  Ci  ■  129  A/V,  V  is  in  volts,  and  I  is  in  amps  [Ref.  14J. 

Current,  pulse  shape,  and  number  of  pulses  varied  from  day  to 
day,  but  consecutive  shots  were  fairly  reproducible. 

2.  Waveform  Capture  by  DCS 

The  radiation  signal  received  by  the  antenna  was  transmitted  to 
the  oscilloscope  by  means  of  the  shielded  coaxial  cable.  The  image  dis¬ 
played  on  the  oscilloscope  was  automatically  captured  by  the  DCS.  The 
signal  captured  by  the  camera  was  immediately  digitized  by  the  frame 
store  board.  The  original  analog  signal  was  displayed  on  a  TV  monitor 
while  the  digitized  signal  was  graphed  on  the  computer  monitor.  Prior  to 
each  Phermex  shot,  the  current  oscilloscope  settings  and  intended 
waveform/file  name  for  each  waveform  were  input  into  the  Compaq  com¬ 
puter  by  use  of  the  DCS  software.  Each  waveform  was  saved  on  the  hard 
disk  of  the  respective  computer  and  a  second  copy  of  each  file  was  pre¬ 
served  on  diskette. 
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m.  DATA  REDUCTION 


A.  DATA  PILE  CONVERSION  AND  TRANSFER 

The  waveforms  were  recorded  and  saved  in  a  binary  format.  Since 
binary  is  the  normal  format  when  viewing  the  computer  display  and 
binary  requires  less  disk  storage,  this  was  the  preferred  format.  As  previ¬ 
ously  mentioned,  various  waveform  manipulations  can  be  accomplished 
using  the  DCS  software  in  the  binary  format.  However,  some  waveform 
manipulations  such  as  the  Fast  Fourier  Transform  (FFT)  cannot  be  per¬ 
formed  with  the  DCS  software.  A  VAX  system  was  selected  to  handle  the 
additional  waveform  manipulations.  Therefore,  an  ASCII  file  format  was 
desirable.  Nine  hundred  and  seventy-five  (975)  binary-formatted  files 
were  individually  converted  to  an  ASCII  format  using  the  DCS  software’s 
built-in  feature.  During  this  process,  each  file  was  coded  with  an  eight¬ 
digit  file  name  to  provide  the  maximum  information  pertaining  to  the  sig¬ 
nal’s  capture  parameters. 

Once  reformatted  and  coded,  all  files  were  transferred  from  the  IBM- 
compatible  computer  to  the  VAX  system.  A  KERMIT  software  program 
provided  the  means  for  transferring  the  files.  The  VAX  system  used  was  a 
Digital  Microvax  computer 

B.  INTERACTIVE  DATA  LANGUAGE  (IDL)  PROGRAM  APPLICATION 

A  program,  written  with  Interactive  Data  Language  (IDL)  [Ref.  19)  on 
the  VAX  system,  was  developed  to  analyze  the  collected  waveforms. 
Appendix  B  contains  a  copy  of  the  main  (modified)  analysis  program. 
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This  program  corrects  for  the  required  in-line  attenuation,  the  antenna 
response  function,  and  the  DC  offset.  Once  the  corrections  are  made,  the 
radiation  signal  and  the  FFT  of  the  signal  are  plotted.  Directly  to  the  right 
of  these  plots,  the  corresponding  B  dot  signal  and  its  FFT  are  depicted. 
Figure  12  is  an  example  of  the  program’s  output.  This  program  was  run 
in  a  batch  mode  to  produce  graphs  for  all  the  radiation  signals.  To  obtain 
additional  information,  the  program  was  modified  to  be  user-interactive. 
This  interactive  program  allows  the  user  to  truncate  the  radiation  or 
B  dot  signal  to  the  exact  start  and  finish  of  the  signal.  In  addition,  the 
program  will  compile  the  peak  frequencies  of  the  FFT  and  their  corre¬ 
sponding  relative  intensities. 

C.  ANTENNA  RESPONSE 

Two  of  the  three  TEM  horns,  A  and  Bl,  were  calibrated  by  the 
National  Bureau  of  Standards  (NBS).  Two  sets  of  responses  for  each  TEM 
horn  were  tabulated.  Each  print-out  of  the  data  set  consisted  of  101 
points.  The  first  set  contained  responses  from  50  MHz  to  5  GHz  at  50 
MHz  intervals  and  a  graph  of  the  data  (Figure  13).  The  second  set  con¬ 
sisted  of  responses  from  1  MHz  to  101  MHz  at  1  MHz  intervals  and  a 
graph  of  the  responses  from  1  MHz  to  512  MHz  (Figure  14).  The  overall 
uncertainty  of  the  calibration  is  ±1.0  dB. 

A  print-out  with  the  response  for  eveiy  MHz  from  1  MHz  to  600  MHz 
would  have  been  ideal.  However,  an  attempt  was  made  using  the  second 
set’s  graph  (Figure  14)  to  extract  the  response  from  5  MHz  to  512  MHz  at 
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RELATIVE  INTENSITY  SIGNAL  (Volts) 


Figure  12.  Example  of  Program  Output 
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A.  1 


Frequency  (MHz) 


Figure  14.  Antenna  Response 
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5.3  MHz  intervals.  This  was  referred  to  as  the  “fine  correction."  The  val¬ 
ues  provided  by  the  calibration  print-out  corresponding  to  Figure  9  were 
referred  to  as  the  “coarse  correction.” 

The  correction  values  are  referred  to  by  the  NBS  National  Engineer¬ 
ing  Laboratory  [Ref.  201  as  the  TEM  horn  transfer  function.  The  transfer 
function  is  defined  as 


TF  =  20  LOG  g 


(8) 


where 

V  =  output  voltage  of  TEM  horn,  rms  volts,  and 
E  =  electric  field  strength,  volts/meter 
The  transfer  function  correction  using  equation  8  was  included  in 
the  analysis  program  to  calculate  the  rFTs.  Points  lying  between  the  5 
MHz  interval  in  the  fine  correction  case  and  between  50  MHz  interval  in 
the  coarse  correction  case  were  derived  by  linear  interpolation. 

A  comparison  of  the  fine  correction  to  the  coarse  correction  is  seen 
in  Figure  15.  The  uncorrected  FFT  is  also  plotted  to  show  the  effect  of  the 
transfer  function.  The  uncorrected  FFT  displays  the  characteristic  50 
MHz  peak  and  its  harmonics.  The  coarse  correction  has  a  similar  trace 
but  it  was  increased  in  intensity  by  a  factor  of  about  100.  However,  the 
fine  correction  causes  a  flattening  of  the  characteristic  peaks.  This  effect 
is  attributed  to  rapid  fluctuations  in  the  graph  of  responses  and  to  the 
method  of  evaluating  the  closest  value.  The  fine  correction  had  an  added 
disadvantage:  it  lengthened  the  program  so  much  that  it  could  not  be 
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executed  in  IDL.  The  coarse  correction  was  therefore  used  in  evaluating 
all  the  FFTs.  Figure  16  depicts  the  transfer  function  for  both  antenna  A 
and  antenna  Bl. 

Also  from  the  NBS  test  procedure,  the  voltage  received  is  related  to 
the  electric  field  E  at  the  antenna,  by 

E  *  CfV  (9) 

where  Cf  =  96.6,  V  is  in  volts,  and  E  is  in  volts/m  (Ref.  201. 
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Figure  16.  Antenna  A  &  El  Transfer  Functions 


■ 
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D.  SIGNAL  REPRODUCIBILITY 

During  the  experiment,  the  reproducibility  of  signals  was  investi¬ 
gated.  A  set  of  three  signals  obtained  at  the  same  angle,  under  the  same 
beam  parameters,  and  with  the  same  antenna  (A)  are  shown  in  Figure 
17.  One  immediately  notices  that  the  three  signals  are  very  similar.  How¬ 
ever.  minor  but  observable  differences  in  the  peaks  are  apparent.  The  dif¬ 
ferences  can  be  attributed  to  slight  changes  in  the  electron  beam  coupled 
with  an  uncalibrated  DCS  system  and  slight  DC  offsets.  In  general,  the 
recorded  signals  for  each  data  set  are  reproducible. 

Figure  18  is  a  set  of  three  signals  taken  under  the  same  conditions 
as  Figure  17.  except  that  antenna  B1  was  used.  Here,  similar  effects  are 
seen.  Again,  the  overall  conclusion  is  that  the  recorded  signals  are  repro¬ 
ducible  in  each  data  set. 

E.  ANTENNA-TO-ANTENNA  COMPARISON 

As  seen  in  Table  1,  antenna  A  was  used  primarily  at  the  smaller 
angles  of  5,  10,  and  15  degrees.  B1  or  B2  was  normally  placed  at  the 
larger  angles  of  30,  45,  and  60  degrees.  Therefore,  a  comparison  of 
antenna  A  to  antenna  B  was  required  to  relate  the  smaller-  and  larger- 
angled  values.  Figure  19  depicts  the  signals  from  antennae  A  and  Bl. 
The  beam  shot  and  angular  position  were  exactly  the  same  for  both 
antennae.  The  signals  display  similar  oscillations  and  characteristics. 
However,  the  maximum  peak-to-peak  voltage  of  antenna  Bl  is  approxi¬ 
mately  one-fifth  of  antenna  A.  Antenna  B2  was  also  compared  to  Antenna 
A  in  a  similar  manner.  Table  2  shows  the  antenna  conversions  for  Bl 
and  B2  to  antenna  A. 
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Figure  19.  Comparison  of  Antenna  A  and  B1  Signals 


TABLE  2 

ANTENNA  CONVERSION 


Antenna 

Factor  to  Convert  Antenna  B 
Signal  to  Antenna  A 

CPC  101  (Antenna  Bl) 

4.9 

CPC  102  (Antenna  A) 

1.0 

CPC  103  (Antenna  B2) 

0.33 

It  was  found  that  antenna  B1  compared  much  more  favorably  than 
antenna  B2  to  antenna  A,  The  signals  from  antenna  B2  were  very  sharp 
and  appeared  almost  erratic.  Antennae  A  and  B1  signals  were  much 
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smoother.  However,  as  mentioned  in  the  previous  section,  the  signals 
were  reproducible  for  each  data  set  recorded  with  the  same  antenna. 
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IV.  RESULTS  AND  DISCUSSION 

.  ■  I  ,  ;  .  ' 

A.  RESULTS 

The  experimental  results  are  organized  into  five  areas:  signals  and 
related  spectrum,  signal  strength  as  a  function  of  angle,  relative  intensity 
of  spectrum  harmonics,  shifting  of  fundamental  frequency,  and  modula¬ 
tion  of  spectrum,  A  discussion  of  the  following  results  appears  in  the 
next  section. 

' 1  \ 

1.  Signals  end  Related  Spectrum 

A  complete  set  of  all  radiation  and  B  dot  signals  along  with  their 
respective  spectra  was  created  using  the  IDL  program.  However,  here  one 
signal  per  angle  (5,  10,  15,  SO,  45.  arid  60  degrees)  and  beam  parameter 
are  provided  for  completeness.  Figure  20-22,  Figures  23-25,  Figures  26- 
28.  and  Figures  29-31  are  signal  and  B  dot  results  for  full  current,  half 
current,  diffuse  beam  (full  current),  and  half  energy  beam  (half  current), 
respectively. 

2.  Signal  Strength  as  a  Function  of  Angie 

For  each  radiation  signal  from  each  data  set,  the  maximum 
peak-to-peak  voltage  was  measured.  This  value  divided  by  two  was  also 
converted  into  the  actual  electric  field  strength  using  equation  9.  The 
current  for  each  beam  shot  was  computed  using  equation  7  and  tabu¬ 
lated  in  Table  3.  Difficulties  with  the  B  dots  on  the  first  day  caused 
inaccurate  current  readings.  It  is  assumed  that  the  current  was 
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Figure  20.  Signal  and  Spectrum  for  Full  Current,  Full  Energy, 
Angles  of  5  and  10  Degrees 
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Figure  22.  Signal  and  Spectrum  for  Full  Current,  Full  Energy, 
Angles  of  45  and  60  Degrees 
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Figure  23.  Signal  and  Spectrum  for  Half  Current,  Full  Energy, 
Angles  of  5  and  10  Degrees 
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Figure  25.  8ignal  and  Spectrum  for  Half  Current,  Full  Energy, 
Angles  of  45  and  60  Degrees 
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Figure  27.  Signal  and  Spectrum  for  Diffuse  Beam,  Full  Current, 
Full  Energy,  Angles  of  15  and  30  Degrees 
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Figure  28.  Signal  and  8pectram  for  Diffuse  Beam,  Full  Current, 
Full  Energy,  Angles  of  45  and  60  Degrees 
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Figure  29.  8ignal  and  Spectrum  for  Half  Energy,  Half  Current, 
Anglea  of  5  and  10  Degree! 
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Figure  31.  Signal  and  Spectrum  for  Half  Energy,  Half  Current, 
Angle*  of  45  and  60  Degrees 
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TABLE  3 

ACTUAL  BEAM  CURRENTS  (IN  AMPS) 


Beam 

Parameter 

Angles  hi  Degrees  | 

10 

10 

IB 

MESH 

45 

*1  j 

Full  Current 

495  ±  13 

286 ± 12* 

275  ±  i4* 

263  ±  5* 

270  sL  2* 

122291 

Half  Current 

302  ±  15 

404  ±2 

297  ±2 

265  ±  5 

310  ±7 

292  ±  13 

Diffuse 

463  ±  22 

449  ±  5 

456  ±5 

463  ±  22 

449  ±5 

456  ±5 

Half  Energy 

312  ±2 

288  ±  13 

2Q8  ±  13 

303  ±  20 

319  ±  1 

303  ±  7 

•Problems  with  B  dot  registered  In  accurate  currents.  Actual  currents  were  approxi¬ 
mately  500  amps. 


approximately  500  A  for  that  day  (for  full-current  tight  beam).  Since  the 
current,  angle,  and  distance  are  ail  known  values,  the  expected  electric 
field  can  be  determined  using  equation  6.  A  comparison  of  expected  to 
actual  electric  fields  is  in  Table  4,  The  average  radiation  peak-to-peak 
value  was  also  plotted  as  a  function  of  angle.  Full  versus  half  current 
(Figure  32),  tight  versus  diffuse  beam  (Figure  33),  and  full  versus  half 
current  (Figure  34)  were  compared.  In  addition,  the  ratio  of  the  upper 
curve  to  the  lower  curve  for  the  plotted  points  was  compared.  The  ratio  of 
curves  for  full  versus  half  current,  tight  versus  diffuse  beam,  and  full 
versus  half  energy  are  exhibited  in  Figures  35,  36,  and  37,  respectively. 

3,  Relative  Intensity  of  Spectrum  Harmonics 

The  fundamental  frequency  and  harmonics  are  easily  observed 
in  the  spectrum.  By  using  the  interactive  version  of  the  IDL  program,  the 
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TABLE  4 

COMPARISON  OF  ACTUAL  TO  EXPECTED  ELECTRIC  FIELD 


Degree* 

8 

10 

18 

30 

48 

_ 

Full 

Current 

A 

Actual 

'xlO3-') 

14.4  1  .7 

12  ±1 

4.9  ±  .6 

2.01.3 

2.121.06 

2.31.1 

Expected 

^lO3--! 
t  tl \) 

9.85  ±.09 

6.4  1.1 

4.21.1 

2.081.02 

1.351.01 

1.001.01 

Half 

Current 

Actual 

t  m) 

9.85  1  .09 

3.3  ±.l 

2.61.1 

n 

1.51.1 

1.32  1 .05 

Expected 

n 

5.04  ±  .02 

2.46 1 .02 

1.27 1 .03 

.82  1 .02 

.56  1 .00 

relative  intensity  for  each  harmonic  J  was  then  plotted  versus  angle.  Thi  s 
was  repeated  for  the  first  four  higher  harmonics  (Figures  38-41)  for  the 
four  different  beam  parameters.  However,  the  full  current  was  carried  out 
to  the  ninth  harmonic  (see  Figure  42). 

4.  Harmonic  Frequency  Shift 

Again,  by  using  the  interactive  version  of  the  IDL  program,  the; 
frequency  of  each  harmonic  was  evaluated.  The  average  harmonic  for 
each  shot  was  calculated  by  plotting  the  frequency  against  the  harmonic 
number  J.  A  least  square  fit  of  the  points  provided  the  slope  or  the 
overall  average  frequency.  The  slopes  for  shots  at  the  same  angle  and 
beam  parameters  were  averaged.  This  average  frequency  (approximately 
50  MHz)  was  then  plotted  versus  angle  for  both  full  (Figure  43)  and  half 
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Figure  32.  Signal  Strength  of  Fuli  Versus  Half  Current  (Full  Energy) 
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Figure  33.  Signal  Strength  of  Tight  Versus  Diffuse  Beam 
(Full  Current,  Full  Energy) 
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Figure  35.  Ratio  of  Full  to  Half  Current  Signals 
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Figure  36.  Ratio  of  Tight  to  Diffuse  Beam  Signals 
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Figure  37.  Ratio  of  Full  to  Half  Energy  Beam  Signals 


51 


RELATIVE  INTENSITY  (x  1000)  RELATIVE  INTENSITY  (x  1000) 


RELATIVE  INTENSITY  (x  1000)  RELATIVE  INTENSITY  (x  1000) 


RELATIVE  INTENSITY  (y  1000) 


-  • .  I'ff  .v,  '  \T  S*'!  1 


r~ . tt— "  . . . 

i 

■  HALc  energy 

I  -  V 


0  ^0  :0  JO  *0  50  80  70 

ANGLE  (Degress) 


HA'_r  energy 
J  «  2 


0  ‘0  20  JO  *0  so  60  70 

ANGLE  (Degrees) 


Figure  40.  First  Four  Harmonies  of  Half  Energy  Spectrum 

(Half  Current) 
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Figure  42.  5  -0  Harmonics  of  Full  Current  Spectrum  (Full  Energy) 
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Figure  43.  Average  Frequency  of  Spectrum  for  Full  Current 

Figure  44)  current.  The  same  calculations  were  performed  on  the  B  dot 
spectrum,  ?nd  the  respective  plots  are  in  Figure  45  and  46.  To  determine 
whether  shifting  of  the  harmonic  frequency  had  occurred,  the  B  dot  full 
current  average  frequency  (Figure  45)  was  subtracted  from  the  full  cur¬ 
rent  average  frequency  (see  Figure  47).  Similarly,  the  difference  was  eval¬ 
uated  and  plotted  for  the  half  current  values  (Figure  48). 

5.  Modulation  of  Spectrum 

The  FFT  results  of  the  B  dot  are  similar  to  the  results  of  the  pat¬ 
tern  seen  from  a  diffraction  grating  in  optics.  From  optics  [Ref.  2IJ,  the 
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Figure  46.  Average  Frequency  of  Spectrum  for  Half  Current  B  Dot 
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Figure  47.  Signal  and  B  Dot  Difference  for  Full  Current 


where  c  Is  the  speed  of  light  and  d  is  the  slit  width  or  pulse  width.  A  clear 
example  of  modulation  is  observed  in  Figure  49.  Here,  a  simulated  B  dot 
pulse  is  formed  by  seven  equal-volt  2.6  ns  pulses  separated  by  20  ns.  In 
this  case,  the  pulse  width  (2.6  ns)  acts  as  the  slit  width  and  is  related  to 
the  modulation  minima.  The  corresponding  minimum  is  therefore  the 
reciprocal  of  2.6  ns,  which  is  365  MHz.  The  20  ns  pulse  separation  forms 
the  50  MHz  peaks  and  higher  harmonics. 
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Figure  49.  Modulation  Example  from  a  Perfectly  Pulsing  Accelerator 


Modulation  minima  were  also  observed  in  the  actual  spectrum  of  the 
signal  and  B  dot.  An  actual  B  dot  and  its  spectrum  are  seen  in  the  top 
right  section  of  Figure  24.  Although  the  beam  is  not  perfect,  it  exhibits 
similar  modulation  minima  at  a  slightly  lower  frequency  (360  MHz).  Table 
5  lists  first,  second,  and  third  modulation  minima  for  various  angles  and 
beam  parameters.  These  values  have  up  to  a  five  percent  margin  of  error. 
In  addition,  only  those  signals  listed  in  Table  5  had  clear  enough  modu¬ 
lation  minima  to  discern  a  value. 
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TABLE  5 

MODULATION  MINIMA 


Beam 

Parameter 

1st 

Minimum 

2nd 

Minimum 

3rd 

Minimum 

Full 

Current 

B  Dot 

275 

525 

750 

5° 

270 

450 

610 

B  Dot 

275 

580 

905 

10° 

280 

— 

— 

B  Dot 

285 

595 

890 

60° 

290 

600 

895 

Half 

Current 

B  Dot 

360 

705 

895 

15° 

320 

675 

— 

B  Dot 

275 

500 

720 

60° 

295 

600 

895 

Half 

Energy 

B  Dot 

220 

470 

740 

5° 

220 

475 

720 

B  Dot 

220 

470 

740 

10° 

240 

450 

650 

B  Dot 

230 

505 

795 

15° 

235 

495 

755 

B  Dot 

— 

445 

860 

30° 

— 

440 

815 

B  Dot 

— 

445 

795 

60° 

— 

450 

740 

B.  DISCUSSION 

The  five  areas  outlined  in  the  results  section  will  now  be  discussed. 
Some  areas  overlap  and  are  inserted  wherever  clarification  is  necessary. 
1.  Signals  and  Related  Spectrum 

There  are  a  number  of  general  observations  which  hold  for  all 
signals  independent  of  antenna  angle  or  beam  parameter.  First,  the 
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overall  pulse  train  length  of  the  radiation  signal  and  the  B  dot  signal  are 
essentially  equal.  Second,  since  the  radiation  depends  on  the  rate  of 
change  of  the  current,  there  Is  a  one-to-one  correspondence  of  B  dot 
,  peaks  to  major  radiation  signal  oscillations.  One  can  also  observe  the 

successive  differences  in  B  dot  pulse  peaks  as  a  similar  difference  in  the 
successive  peak  oscillations.  Last,  the  50  MHz  frequency  and  higher 
harmonics  are  clearly  evident  in  all  radiation  spectra.  However,  the  dif¬ 
fuse  beam  radiation  spectrum  harmonics  were  less  consistent  and  dis¬ 
tinguishable.  Also,  the  modulation  minima  of  the  diffuse  beam  were  not 
clearly  evident  from  the  spectrum. 

2.  Signal  Strength  at  a  Function  of  Angle 

By  viewing  Figures  32,  33,  and  34,  one  discovers  a  general 
decrease  in  signal  strength  with  increase  in  angle.  This  effect  is  mainly 
attributed  to  transition  radiation.  A  comparison  of  the  general  transition 

*  radiation  pattern  (Figure  3)  and  the  radiation  (Figures  32-34)  observed 
from  this  experiment  shows  significant  similarities.  Actual  electric  field 
strength  from  Table  4  is  comparable  to  the  expected  strength  at  angles  of 
15  degrees  or  greater.  However,  at  smaller  angles,  the  actual  electric  field 
is  considerably  larger.  This  increase  in  field  strength  may  be  attributed  to 
diffracted  Cherenkov  radiation  and  diffracted  transition  radiation.  A 
closer  look  at  the  5-  and  10-degree  values  of  both  curves  on  Figure  34 
shows  a  large  difference  at  these  angles.  The  full  energy  (29.5  MeV)  is 
above  the  energy  to  cause  Cherenkov  radiation,  and  the  half  energy  (15 
MeV)  is  below  the  required  energy.  The  ratio  from  Figure  37  at  small 

•  angles  also  suggests  diffracted  Cherenkov  radiation  may  be  present. 
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The  ratio  of  the  full  current  peak-to-peak  signal  versus  the  half 
current  peak-to-peak  signal  (Figure  35)  exhibits  (with  the  exception  of  the 
10-degree  data  point)  a  fairly  constant  ratio  of  approximately  1.4.  On  re¬ 
examining  equation  6.  one  finds  that  the  ratio  of  full  versus  half  current 
signals  at  the  same  angle  is  merely  the  ratio  of  the  currents.  All  other 
values  in  equation  6  remain  constant  and,  therefore,  from  Table  5  the 
average  ratio  of  full  to  half  currents  is  approximately  1.6. 

In  Figure  33,  a  propagating  (tight)  beam  and  a  non-propagating 
(diffuse)  beam  are  compared.  Again,  one  observes  a  similar  transition 
radiation  pattern.  I-ooking  again  at  Figure  36,  one  notices  a  small  ratio 
(of  approximately  2)  for  the  smaller  angles  and  a  larger  ratio  (of  approxi¬ 
mately  8)  for  the  larger  angles.  It  appears  that  a  non-propagating  beam 
still  has  significant  radiation  in  the  forward  direction  but  radiates  signifi¬ 
cantly  less  than  a  propagating  beam  in  the  other  directions. 

3.  Relative  Intensity  of  Spectrum  Harmonics 

As  previously  mentioned,  the  50  MHz  harmonics  are  readily 
observable  in  the  radiation  spectrum.  Figure  38-41  Indicate  the  relative 
intensity  of  50  MHz  (J  ■  1)  and  the  next  three  higher  harmonics.  In  gen¬ 
eral,  the  harmonic  intensity  follows  the  same  pattern  as  the  signal 
strength.  Therefore,  increases  in  angle  cause  a  decrease  in  relative 
intensity.  Also,  each  successive  harmonic  curve  decreases  in  intensity. 
The  relative  intensity  is  still  significant  in  higher  order  (sometimes  up  to 
J  =  20)  harmonics.  This  effect  is  seen  in  each  signal  FFT  and  is  graphi¬ 
cally  depicted  in  Figure  42  with  J  *  5  through  9. 
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4.  Harmonic  Frequency  Shift 

Figures  45  and  46  indicate  that  a  relatively  constant  50  MHz 
beam  was  being  produced.  Viewing  Figures  43  and  44,  a  frequency  shift 
with  increasing  angle  is  indicated.  To  negate  any  effect  changes  in  the 
beam  source  had  on  this  shifting,  the  average  B  dot  frequency  was  sub¬ 
tracted  from  the  average  signal  spectrum  frequency  (Figures  46  and  47). 
Again,  a  slight  shift  to  higher  frequency  with  increasing  angle  is  evident. 
This  predicted  pattern  was  seen  earlier  in  Figure  2.  These  experimental 
findings  seem  to  support  the  theory  behind  Figure  2.  However,  measure¬ 
ments  in  both  Figures  43  and  44  involve  two  different  collection  systems. 
One  collection  system  captured  the  signals  at  5,  10,  and  15  degrees, 
whereas  the  other  system  collected  the  signals  at  the  remaining  angles. 

5.  Modulation  of  Spectrum 

The  summary  of  modulated  minima  (Table  5)  excludes  minima 
from  the  diffuse  beam.  As  mentioned  previously,  modulation  minima 
were  difficult  to  locate  for  diffuse  beams  and  were  consequently  omitted. 
At  larger  angles,  the  modulation  minima  were  easier  to  discern.  This  fact 
is  in  agreement  with  the  BuskirK,  Neighbours,  and  Maruyama  model 
(Ref.  91.  From  the  remaining  beam  parameters,  there  is  significant 
agreement  between  the  B  dot  and  signal  spectrum  minima.  In  general, 
the  first  and  sometimes  the  second  minima  are  consistent.  The  first  min¬ 
ima  appears  on  the  average  at  235  MHz  for  a  five-pulse  beam,  280  MHz 
for  a  six-pulse  beam,  and  340  MHz  for  a  seven-pulse  beam.  Using  equa¬ 
tion  10  and  multiplying  by  the  respective  number  of  beam  pulses 
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consistently  gives  values  between  3.4  and  3.6  ns.  These  values  are  very 
close  to  the  actual  value  of  the  beam  pulse  width  (3.5  ns). 
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V.  CONCLUSIONS 


Radio  frequency  radiation,  produced  by  the  electron  beam  of  the 
PHERMEX  accelerator,  consisted  mainly  of  transition  radiation  (EMP) 
and  minor  contributions  from  diffracted  Cherenkov  radiation  (when  at 
full  energy).  Regardless  of  the  beam  energy,  current,  or  frequency,  the 
following  general  observations  can  be  made  concerning  the  analyzed 
radiation  and  its  associated  spectrum: 

1.  There  is  a  one-to-one  relation  between  the  radiation  oscillating 
peaks  and  the  electron  beam  pulses. 

2.  The  length  of  the  radiation  pulse  train  remained  constant  and  con¬ 
sistent  with  the  beam  pulse  train. 

3.  PHERMEX  operating  frequency  of  50  MHz  was  consistently 
observed  with  up  to  the  20th  harmonic  in  the  signal  spectrum. 

4.  In  addition,  there  is  some  evidence  that  shifting  may  occur  to  the 
fundamental  frequency  and  its  harmonics  in  the  signal  spectrum. 
Increasing  angle  seems  to  cause  an  increase  in  the  average  har¬ 
monic. 

5.  Electric  field  strength  of  the  radiation  is  proportional  to  the  beam 
current. 

6.  Modulation  minima  of  radiation  signal  spectrum  were  similar  to  B 
dot  spectrum.  The  first  modulation  minima  is  directly  related  to  the 
beam  pulse  width. 

In  addition,  differences  between  propagating  and  non-propagating 
beams  were  noted  as  follows: 

1.  Propagating  beams  retain  significant  beam  characteristics  through 
their  spectrum,  whereas  non-propagating  beams  do  not. 

2.  Non-propagating  beams  produce  less  radiation  in  the  forward  direc¬ 
tion  than  a  propagating  beam. 
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3.  A  significant  decrease  in  radiation  exists  between  propagating  and 
non-propagating  beams  in  other  than  the  far  forward  direction. 

There  has  been  a  major  improvement  in  the  collection  system.  Sig¬ 
nals  can  now  be  quickly  digitized  and  stored  in  a  computer  for  later  anal¬ 
ysis.  Previous  experiments  resulted  in  only  analog  measurements. 

Further  work,  both  theoretical  and  experimental,  is  already  planned 
to  continue  investigating  this  and  additional  data  collected  during  the 
PHERMEX  experiment. 
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APPENDIX  A 

EQUIPMENT  LIST 


1.  Tektronix  7103  Oscilloscope  (3  each) 

2.  Transverse  Electromagnetic  Horn  Antenna  (3  each) 

3.  Tektronix  Digitizing  Camera  System  (DCS)  (3  each) 

4.  Personal  Computers 

Compaq  Portable  II  Model  4  (2  each) 

IBM  XT  (1  each) 

5.  Video  Monitors  (3  each) 

6.  Double- Shielded  Cables  (3  each) 

7.  Tektronix  Video  Processor  HCOZ 

8.  IBM  Graphics  Printer 
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APPENDIX  B 

IDL  ANALYSIS  PROGRAM 


;  This  is  an  updated  version  of  rf_sxpt  by  Rich  Lally  on  9  MR  90. 

; PROGRAM  TO  ANALYZE  RF  DATA  TAKEN  WITH  DCS 
FOR  X-1,S  DO  BEGIN 
PRINT, '  ' 

CNDFOR 

PRINT,”  PROGRAM  PLOTTER. PRO"  /Program  Titls 

PRINT,"  " 

PRINT, “***•  Program  to  analyse  DCS  Data.  ***** 

FOR  1-1,5  DO  BEGIN  ; Skips  5  lines 

PRINT, '  ' 

ENDFOR 

READ, 'DO  YOU  NISH  TO  FIND  MAXIMUMS  ?  1  FOR  YES  AND  2  FOR  NO',FMAX 


PRINT, "File  MOST  NOT  have  an  .ASC  extension." 

;  .mmmmmmmmmmm  INPUT  OF  FILENAME  mmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmmm 

NAME-' ' 

READ, 'Enter  the  name  of  the  data  file', NAME  .Resets  file  name  to  null 

/  Set 

FOR  JJ-1,2  DO  BEGIN 


/Measures  length  of  file  name 


; INPUTS  CORESPONDING  BDT  NAME 
; Stripa  off  direatory  t  . asc 

DAYLETTER-STRMID  (TITLE,  0,1) 

CASE  DAYLETTER  OF  * 

'M' !  DAY-' MON' 

'T':  DAY**'  TUE’ 

'H' :  DAY-' NED' 

'H'l  DAY-' THU' 

'F'  :  DAY-' mi' 

ELSE :  DAY-' UNXNONN ' 

ENDCA8E 


NAMELEN— STRLEN  (NAME) 

IF  JJ  IQ  2  THEM  8TRVUT, NAME, 'BOX', 1 
TITLE— STRMID  (NAME,  0,  (NAMELEN)  ) 


DEGREES— STRMID (TITLE, 1,2) 
ANTENNA— STRMID (TITLE, 3,  1) 
SMOTNUMBER— STRMID (TITLE, 4,3) 
CURRENTL-STRMID (TITLE , 7 , 1 ) 

CASE  CURRXMTL  OF 

'F' :  CURRENT-' FULL' 

l  CURRENT-' HALF ' 

ELSE:  CURRENT-' UNKNOWN' 

SNDCASE 


IF  DEGRESS  EQ  'BD'  TEEN  TOPLINE-TITLE  ♦  '-'  *  CURRENT  ♦  '!,'» 
*  'B  DOT  FOR  *  +DAY  SMOTNUMBER  ELSE  $ 

TOPLINE-TITLE  ♦  ♦  CURRENT  +  '  I,’  *  DEGREES  *  '  DEG, ANT  '  * 

*  ANTENNA  ♦  ','  +OAY  ♦  '  #'  ♦  SMOTNUMBER 

PILE  -  NAME  *  ' .ASC' 

■  — — —  READ  DATA  IN  PILE  ■  — . .  .i-.— ■ 
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; This  is  an  updated  version  of  rf_expt  by  Rich  bally  on  8  APR  90. 

; PROGRAM  TO  ANALYZE  RF  DATA  TAKEN  WITH  DCS 
TOR  1-1,5  DO  BEGIN 
PRINT, '  ' 

ENDrOR 

PRINT,"  PROGRAM  PLOTTER. PRO"  /Program  Title 

PRINT, "  " 

PRINT,"****  Program  to  analyze  DCS  Data.  ****■• 

FOR  1-1,5  DO  BEGIN  ; Skips  5  linas 

PRINT, '  ' 

ENDFOR 

READ, 'DO  YOU  WISH  TO  FIND  MAXIMUMS  ?  1  FOR  YES  AND  2  FOR  NO',FMAX 


PRINT, "Fila  MUST  NOT  hava  an  .ASC  extension." 

;  mmmmmmmmmmmm  INPUT  OF  FILENAME  ————————————— 

NAME-' ' 

READ, ' Entar  tha  nama  of  tha  data  fila', NAME  ;Raaats  fila  nama  to  null 

/  Sat 

FOR  JJ-1,2  DO  BEGIN 


/Maaauras  length  of  fila  nama 


/INPUTS  CORESPONDING  BDT  NAME 
/Strips  off  diractory  &  . asc 

D-.iLETTER-STRMID  (TITLE,  0,  1) 

CASE  DAYLETTER  OF 
' M' :  DAY-' MON' 

'T':  DAY—' TUE' 

'W' :  DAY-' WED' 

'H':  DAY-' THU' 

'F'!  DAY^'FRI' 

ELSE :  DAY- ' UNKNOWN ' 

ENDORSE 


NAMELEN-STRLEN (NAME ) 

IF  JJ  EQ  2  THEN  STRPUT, NAME, ' BDT' , 1 
TITLE-STRMIO (NAME, 0, (NAMELEN) ) 


DEGREES-STRMIU (TITLE, 1,2) 

ANTENNA— STRMID (TITLE -3,1) 

SHOTNUMBER-3TRMID (TITLE, 4, 3) 

CURRENTL-STRMID (TITLE ,7,1) 

CASE  CURRSNTL  OF 
'F':  CURRENT-'  FULL' 

'H':  CURRENT-' HALF' 

ELSE :  CURRENT- ' UNKNOWN ' 

END CASE 


IF  DEGREES  EQ  'BD'  THEN  TOPLINE-TITLE  +  '-'  +  CURRENT  +  '  I,’  S 
+  'B  DOT  FOR  '+DAY  +  '  #'  +  SHOTNUMBER  ELSE  $ 

TOPLINE-TITLE  4-  '-'  4-  CURRENT  +  '  I,'  ♦  DEGREES  '  DBG,  ANT  '  $ 
+  ANTENNA  *  ' ,  '  4-DAY  4-  '  #'  4-  SHOTNUMBER 

FILE  -  NAME  4-  '  .ASC' 
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OPENR, 2 , FILE  ; Opana  input  til a  for  raiding 

haadar-strarr <80, 30)  ; Oinanaiona  string  arrays 

r dhdr-atrar r (80,1) 

DATA  -  FLTARR(512) 

AA  «  '  ' 

;Raad  in  valuas  from  haadar 


1-0 

WHILE  NOT  EOF  (2)  DO  BEGIN 

IF  (<I  GT  19)  AND  (I  LT  532))  THEN  BEGIN 
READF ,2, '  $ (F8.6) ' ,  FF 
DATA  (1-20)  -  FF 
END IF  ELSE  BEGIN 
READF, 2,  AA 
HEADER ( I )  -  AA 
ENDELSE 
I-I  +  1 
END WHILE 


' :')  +2 

(HEADER (6)  , 

POS1, 

3)  , 

2) 

MM  +2 

(HEADER (7)  , 

POS2, 

3), 

2) 

POS1  -  STRFOS  (HEADER (6) , 

STARLBL  -  STRTRIM  (STRMID 
POS2  -  STRFOS  (HEADER (7), 

STOPLBL  -  STRTRIM  (STRMID 

START  -  FIX (STARLBL)  6  STOP  «  FIX (STOPLBL) 
DELTA-STOP -START  +  1 
POS3  -  STRFOS  (HEADER (12) ,  ':')  +2 
TXMLBL  «  STRTRIM  (STRMID (HEADER (12) ,POS3, 20) ,2) 
TIME  -  FLOAT (TIMLBL) 

IF  (DEGREES  EQ  '  BD' )  AMD  (TIME  G1‘  3.E-9)  THEM  TIME 


; STRING  TO  NUMBER  CONVERSION 


4.34783E-10 


T-indgan (dalta) *tima  /TIME  DATA 

;  ——.—————ATTENUATION  CORRECTION— 

DB-0 

IF  (DAYLETTER  EQ  'MM  AMD  (AMTKMMA  EQ  'A')  THEM  DB-100.0 

IF  (DAYLETTER  NE  'MM  AMD  (ANTENNA  EQ  'AM  THEN  DB-63.1 

IF  (DAYLETTER  EQ  'MM  AMD  (ANTENNA  EQ  'BM  THEN  DB-3.2 

IF  (DAYLETTER  EQ  'TM  AND  (ANTENNA  EQ  'B'J  THEN  DB-1C.0 

IF  (DAYLETTER  EQ  'W')  AMD  (ANTENNA  EQ  'B')  THEN  DB-10.0 

IF  (DAYLETTER  EQ  'F')  AND  (ANTENNA  EQ  'BM  THEN  DB-316.2 

IF  (DAYLETTER  EQ  'H')  AND  (ANTENNA  EQ  'B')  THEN  BEGIN 

IF  (SHOTNUMBER  LT  019)  THEN  DB-10.0 

IF  (f'HOTNUMBCR  GE  019)  AND  (SHOTNUMBER  LT  027)  THEN  DB-100.0 
IF  ( SHOTNUMBER  GE  027)  THEN  DB-316.2 
END IF  ELSE  BEGIN 

IF  (DEGRESS  EQ  'BD')  THEN  DB-70.711 

ENDELSE 


cloaa, 2 


Ixtitla-' TIME  (naae) ' 

! ytitla-' SIGNAL  (Volts) ' 

! type-  16 
I MTITLE-TOFLXNE 

X-T»l . 0E9  /SCALE  TO  NANOSECONDS 

UY-DATA (START : STOP ) 

Y-UY  *  DB  /ATTENUATION  CORRECTION 

;  ———CURSOR  ROUTINE  TO  SIZE  WAVEFORM— — — — — ■ 

SET  PLOT, 1 

IF  Jj  EQ  2  THEN  GOTO, LIFE 
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!C-0 


PLOT, X, T 

GOTO, LIFER 

LIFE:  LX-X  &  LF-F 

ERASE 

! Jctitl*-' TIKE  <na«c) ' 
lytitla-' SIGNAL  (Volta)' 
!typ«-  16 
I MTITLE-TOPLINE 
!XMAX-MAX(X) 

!C-0 
! GRID-0 
IXTICKS-0 


PLOT, LX, LF 

LIFER: EFEEEEE-0 

CURSOR, SI, S2, XI, FI 

CURSOR, S3, S4,X2,F2 

CORRECT-TIME *1 . 0E09 

X2-X2 /CORRECT 

Xl-Xl /CORRECT 

IF  XI  LE  0  THEN  XI  -  0.0 

XXX-X (XI : X2 )  S  FFF-F (XI : X2 ) 

DCOFFSET— F (XI) 

F-FFF-DCOFFSET 
ZEROER-X (XI) 

X-XXX-ZEROER 
DELTA-N_ELEMENTS (X) -1 

;  —————————  FOURIER  TRANSFORM 

I-INDGEN (DELTA) 

NU-1/ (TIME-DELTA)  ; FUNDAMENTAL  FREQUENCF 

P2-DELTA/2 

FREQ-FLTARR (P2 ) 

CDB»FLTARR(P2) 


FOR  J-l, P2-1  DO  BEGIN 
FREQ (J) -I ( J) *NU/1 . 0E6 
BAD-WHERE  (  FREQ  LE  600  ) 
BAD-WHERE  (  FREQ  LE  1200  ) 


;  IN  MHZ 

;TO  CHANGE  MAX  FREQ  CHANGE  NUMBER  I.E.  600 
»•  TO  CHANGE  MAX  FREQ  CHANGE  NUMBER  I.E.  600 


'ANTENNA  CORRECTION' 


IF  DEGREES  EQ  'BD'  THEN  GOTO, DIVE 
IF  ANTENNA  EQ  'A' THEN  BEGIN 


CASE  1  OF 

(FREQ ( J)  GE  1)  AND  (FREQ ( J)  LT- 50) :  CDB(J)»< (FXEQ(J) -1) / (50-  $ 

1) >*(30.1-30.0) +30.0 

(FREQ (J)  GE  50)  AND  (FREQ ( J)  LT  100):  CDS ' J) - ( (FREQ ( J) -50) / (100-  $ 
50))*(30. 3-30.1) +30,1 

(FREQ ( J)  GE  100)  AND  (FREQ (J)  LT  150):  CDB(J)-((FREQ(J)-100)/(150- 
100>) *<37.9-30.3>+30.3 

(FREQ ( J)  GE  150)  AND  (FREQ ( J)  LT  200):  CDB (J) «( (FREQ (CO -150) /  (200- 
150)) *(41.9-37.9! +37 , 9 


(FREQ ( J)  GE  200)  AND  (FREQ (7) 
200) ) * (32 .1-41 . 9) +41 . 9 
(FREQ ( J)  GE  250)  AND  (FREQ (J) 
300) ) * (30 . 7-32 . 1) +32 . 1 
(FREQ (0 )  GE  300)  AND  (FREQ ( J) 
300) ) *(36.0-30.7)+30.7 
(FREQ ( J)  GE  350)  AND  (FREQ(J) 
350) ) * (34. 8-36.0) +36.0 
(FREQ (J)  GE  400)  AND  (FR1Q(J) 


LT  250):  CDB (J) -( (FREQ (J) -200) / (250- 
LT  300):  CDB <J>-( (FREQ (J) -250)/ <300- 
LT  350):  CDB (J)-( (FREQ (J) -300)/ (350- 
LT  400):  CDB (J)-( (FREQ (J) -350)/ (400- 
LT  450):  CDB (J) -( (FREQ (J) -400) / (450- 


$ 

S 

S 

$ 

$ 

$ 

$ 
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400) ) *  <37 . 3-34 .8) +34 . 8 

(FREQ(J)  GE  430)  AND  (FREQ  (J)  LT  300):  CDB  ( J)  -  ( (FREQ(CT)  -450)  /  (300-  $ 

430)  )  * (3S.3-37.3)+37.3 

(FREQ(J)  GE  300)  AND  (FREQ  (J)  LT  330):  CDB (J) -< (FREQ (J) -500) / <350-  $ 

500) )  *  (42.3-38.5)4-38.5 

(FREQ (J)  GE  530)  AND  (FREQ ( J)  LT  600):  CDB (J) -( (FREQ (J) -530) / (600-  $ 

350) ) * (39 . 6-42.3) +42.5 


ELSE:  CDB (J) -40.0 
END CASE 

END IF  ELSE  BEGIN 

IF  (DAY  EQ  'H' )  OR  (DAY  EQ  'F')  THEN  GOTO, DIFFER 
IF  (DAY  EQ  'W')  AND  (SHOTNOMBER  GE  008)  THEN  GOTO, DIFFER 


CASE  1  OF 

(FREQ ( J)  GE  1)  AND  (FREQ ( J)  LT  50):  CDB (J) « < (FREQ (J) -1) / (50-  S 
1)  )*<30.7-33.3)+33.3 

(FREQ  < J)  GE  50)  AND  (FREQ ( J)  LT  100):  CDB ( J) - ( (FREQ (J) -50) / (100-  $ 
50)) *<33. 4-30. 7) +30. 7 


(FREQ ( J)  GE  100)  AND  (FREQ ( J)  LT 

100) ) * (54.0-33.4) +33.4 

(FREQ ( J)  GE  150)  AND  (FREQ ( J)  LT 

15 0))*(34. 6-54.0) +34 . 0 

(FREQ ( J)  GE  200)  AND  (FREQ ( J)  LT 

200))*(32.8-34.6)+34.6 

(FREQ ( J)  GE  250)  AND  (FREQ ( J)  LT 

300) )  * (35 . 0-32 .8) +32 . 8 

(FREQ ( J)  GE  300)  AND  (FREQ ( J)  LT 

300) ) * (37.5-35.0) +35.0 

(FREQ (J)  GE  350)  AND  (FREQ (J)  LT 

330)) *(33.3-37.5) +37.5 

(FREQ ( J)  GE  400)  AND  (FREQ(J)  LT 

400) )  * (37 .4-33.3) +35.3 

(FREQ ( J)  GE  450)  AND  (FREQ (J)  LT 

450) ) * (38 . 4-37 . 4) +37 . 4 

(FREQ ( J)  GE  500)  AND  (FREQ ( J)  LT 

500))*<42.3-38.4>+38.4 

(FREQ (J)  GE  550)  AND  (FREQ ( J)  LT 

550) ) * (41. 7-42. 3) +42. 3 


150):  CDB (J)«< (FREQ (J) -100)/ (150- 
200):  CDB <J)-( (FREQ (J) -150)/ (200- 
250)  :  CDB  (.7)  ■  <  (FREQ  ( J)  -200)  /  (250- 
300) :  CDB ( J) - ( (FREQ ( J) -250) ! (300- 
350):  CDB (J)-< (FREQ (J) -300)/ (350- 
400):  CDB (J) -( (FREQ (J) -330)/ (400- 
450):  CDB <J)-< (FREQ (J) -400)/ (450- 
500):  CDB (J)-( (FREQ (J) -450)/ (300- 
350):  CDB (J) -< (FREQ (J) -500)/ (550- 
600):  CDB (J)«< (FREQ (J) -550)/ <600- 


$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 


ELSE:  CDB  < J) *40 . 0 
ENDCASE 


ENDELSE 
GOTO, FLIP 
DIPPFR:  CASE  1  OF 

(FREQ ( J)  GE  1)  AND  (FREQ(J)  LT  50):  CDB ( J) - < (FREQ ( J) -1) / (50-  $ 

1) )*(51.0-31.5)+31.5 

(FREQ ( J)  GE  50)  AND  (FREQ (J)  LT  90):  CDB ( J) * ( (FREQ ( J) -50) / (90-  $ 

50) ) * (44.5-51 . 0) +51 . 0 

(FREQ ( J)  GE  90)  AND  (FREQ (J)  LT  100):  CDB (J) -( (FREQ (J) -90) / (100-  $ 

90)) *(48. 0-44. 3) +44. 5 

(FREQ (J)  GE  100)  AND  (FREQ ( J)  LT  200):  CDB (J) -( (FREQ (J) -loO) / (200-  $ 

100) >*(53.5-48.0! +48.0 

(FREQ (J)  GE  200)  AND  (FREQ (J)  LT  300):  CDB (J) •( (FREQ (J) -200) / (300-  $ 

200) )  * (33.7-53.5) +55.5 

(FREQ ( J)  GE  300)  AND  (FREQ ( J)  LT  400):  CDB (J) -( (FREQ (J) -300) / (400-  S 
300) )  * (52 . 0-55 .7) +55 .7 

(FRAQ(J)  GE  400)  AND  (FREQ(J)  LT  300):  CDB(O) ■(  (FREQ (J) -400) / (500-  $ 

<00! )  * (47 . 0-52 .0) +32 . 0 

(FREQ(J)  GE  500)  AND  (FREQ ( J)  LT  600):  CDB(J) -( (FREQ (J) -500) / (600-  $ 

500) ) * (41 .5-47 . 0) +47 . 0 

(FREQ  ( J)  GE  600)  AND  (FRJEQ  ( J)  LT  700):  CDB (J) *< (FREQ (J) -600) / (700-  $ 

600) )  * (48.0-41.5) +41.5 

(FREQ (J)  GE  700)  AND  (FRSQ(J)  LT  800):  CDB (J) • ( (FREQ (JJ -700) / (800-  $ 
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700) ) *  (49. 0-48.0) +4 8.0 

(FREQ(J)  GE  800)  AND  (FREQ(J)  LT  900):  CDS ( J) - ( (FREQ (J) -800) / (900-  $ 

800) ) *  (43 .3-49.0) +49.0 

(FREQ  ( J)  GS  900)  AND  (FRXQ(.T)  LT  1000):  CDB (J) -< (FREQ (J) -900) / <1000-  S 
900) ) « (43 . 3-43 .3) +43 . 5 

ELSE:  CDB ( J) ■43.5 


END  CAS  F 

GOTO, FLIP 
DIVE:  COB (J) -0.0 
FLIP:  FFF-1 .  0 
SNDFOR 


FREQ (0) -0.0 

IF  ANTENNA  EQ  'A' THEN  COS (0) -30.0 
IF  ANTENNA  EQ  'S' THEN  CDB  (0-33.3 
BAD IE-MAX (BAD) 

; WITHOUT  ANTENNA  RESPONSE 
; ANTENNA  RESPONSE  FUNCTION 


; Z-FFT  (1,1) 
CT-10*  (CDB/20) 


UZ-FFT  (Y,  1) 
Z-UZ*CT 


/WITH  ANTENNA  RESPONSE 
; WITH  ANTENNA  RESPONSE 


MAKE  REGIS  PLOT  OF  DATA  — 1 


SET  PLOT, 3 

IF  Jj  EQ  2  THEN  (NOERAS-1 
1PSYM-0 

IF  JJ  EQ  1  THEN  OPENW, 5, TITLE  +'  COMPOSITE. PLT/Non«' 
PLOT  TO, -5 

IF  EQ  1  THEN  SET  VIEWPORT, . 08, . 48, . 60, . 95 
IF  JJ  EQ  2  THEN  SET~VINWPORT ,  .58,  .98,  .60,  .95 
! GRID-0 

SCALE RR  -  MAX (X) 

IXMAX  -  SCALE RR 
!C-0 

PLOT , X, X 
STOP 


MAKE  REGIS  PLOT  OF  FOURIER  TRANSFORM  — > 
! NOERAS-1 


IF  JJ  EQ  1  THEN  SET  VIEWPORT, .08, .48, .10, .45 
IF  JJ  EQ  2  THEN  SETJ/IEWPORT,  .58,  .98,  .10,  .45 

I  GRID- 1 


!MTITLE-'FFT  OF  '  +  TITLE 

1XTITLE-' FREQUENCY  (MH*) ' 
iYTITLE-' RELATIVE  INTENSITY' 

SCCALER  -  600  /TO  CHANGE  MAX  FREQ  CHANGE  600  TO 

;  WHATEVER  FREQ 

SCCALER  -  1200 
iXKAX-  SCCALER 

SPOT  -  SCCALER/ 3 
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VA1UE-ABS (Z (0iP2-l) ) 


!C-0 

! IGNORE*! 

PLOT  IO,FREQ(0:8ADIE) , VALUE 
; PLOT, FREQ ( 0 : BAD IE ) .VALUE 
i IT  JJ  EQ  X  THEN  STOP 

IF  FMAX  NE  1  THEN  GOTO, FLAGGER 
IF  JJ  EQ  1  THEN  BEGIN 

;  —— —————— DETERMINE 

READ, 'USE  MOUSE  TO  FIND  MAX  VALUES,  BUT  FIRST  ENTER  NUMBER  OF  MAXZMUMS ' , NMAX 
OPENW, 7 , TITLE  +  '  MAXX. DAT/LIST' 

OPENN, 6, TITLE  +  '  ""MAX. DAT/LIST' 

PRINTF , 7 , " 9 <rDATA  VALUES  OF  ' , AS ) " , TITLE 
PRINTF, «,"$ ('DATA  VALUES  OF  ',A8)H, TITLE 
PRINTF, 7,"#  INTENSITY  FREQUENCY" 

PRINTF, 7," 

PRINTF ,  6,  "T (F8  .3Tlr,"NHJK  —  - 

CURSOR,  21, 22,  XX,  Y1 
IF  XI  LE  0  THEN  XX  -  0.0 

FOR  IJ  -  0 , NMAX- 1  DO  BEGIN 
WAIT, 0.5 
X2-X1 

CURSOR, 21, Z2, XI, Y1 
CORRECTO-NU 
X1-X1*1 . OE6/CORRECTO 
IVAL-MAX (VALUE (X2 : XX ) ) 

PRINTF, 7, '$(12, 013.3,013. 5)', I J, IVAL, FREQ (WHERE (VALUE  EQ  IVAL) ) 
PRINTF, 6, I J 
PRINTF, 6, IVAL 

PRINTF , 6 , FREQ (WHERE (VALUE  EQ  IVAL) ) 

ENDFOR 

CLOSE, 6 
CLOSE, 7 

END IF  ELSE  BEGIN 


; PLOTS  SEMI -LOG 
/PLOTS  A  LINEAR  PLOT 


;  — —  mmmmmmmmmmmm -PRINT  MAX  I  MUMS  TO  A  FILE— 

READ, 'USE  MOUSE  TO  FIND  MAX  VALUES,  BUT  FIRST  ENTER  NUMBER  OF  MAXIMUMS '  ,  NMAX 
OPENW, 4, TITLE  +  '  MAXX. DAT/LIST' 

OPENW, 3, TITLE  +  r  MAX. DAT/LIST' 

PRINTF, 4, "$('DXTA  VALUES  OF  ', AS TITLE 
PRINTF, 3, "$(' DATA  VALUES  OF  ', AS) ",  TITLE 
PRINTF, 4,"#  INTENSITY  FREQUENCY" 

PRINTF,  4,"  _  _ " 

PRINTF, 3,  "? (Fa'.TPVNWJOi  ‘ 

CURSOR,  21,  Z2,  XX,  FI 
IF  XX  LE  0  THEN  XX  -  0.0 

FOR  ZJ  -  O.NMAX-l  DO  BEGIN 
WAIT, 0.5 
X2-XX 

CURSOR, Z1,Z2, XX, YX 

CORRECTO-NU 

Xl-XX  * 1 . 0E« /CORRECTO 

IVAL-MAX (VALUE (X2 : XX ) ) 

PRINTF, 4, ' $ (12,6X3. 5, GX3. 5) ' , I J-l, IVAL, FREQ (WHERE (VALUE  EQ  IVAL)) 
PRINTF, 3, IJ 
PRINTF,  3,  IVAL 
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PRINTF , 3 , FREQ (WHERE (VALUE  EQ  IVAL) ) 
ENDFOR 
CLOSE , 3 
CLOSE r  4 

ENDELSS 

FLAGGER :  TRRT-0 
ENDFOR 

AA-.03  *  MIN (VALUE) 

BB“ . 90  *  SCCALER 

XTOUTS , BB, AA, STRMID ( ! STIME, 0 , 11 ) 

CLOSE, 5 
STOP 

CLEAN? LOT 
ERASE 

PLOT_TO,  0 

END 
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